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Tha  vlacoua  supersonic  flew  ovar  a  sharp  cona 
at  lacldaaca  la  examined  numerically  with  a  couplad 
strongly  lapllclt  algorithn  for  tha  propartlaa  In 
tha  plana  normal  to  tha  cona  ax la.  Tha  Nsvier- 
Stokaa  squat Ions  ara  considered  la  a  boundary 
layer-lika  or  parabolized  aaimer  and  global 
relaxation  la  coneldered  for  tha  pressure  Inter¬ 
action.  It  la  shown  that  dapartura  affects  can  be 
effectively  eliminated  by  forward  differencing  for 
the  axial  pressure  gradient.  Moreover,  this 
approximation  retains  the  Implicit  free  pressure 
interaction  required  for  geometries  where  axial 
flow  separation  la  possible. 

A 

1.  Introduction 

Por  a  significant  class  of  geometric  configu¬ 
rations  and  flow  conditions  solutions  of  the  full 
Revier-Stokes  equations  can  be  approximated  quite 
accurately  with  boundary  layer  like  marching 
techniques  in  thin  viscous  layers,  couplad  with 
relaxation  or  marching  procedures  for  lnviacld 
sub-  and  supersonic  regions,  respectively.  Por 
large  Reynolds  number  flows  this  approach  has  been 
developed  la  several  ways,  s.g. ,  Interacting 
boundary  layer  theory,1  parabolised  Havier-Stokas 
(PRS)  theory,2  viscous  or  single  layer  theory3  sad 
two-layer  Interactive  theory.*  Single  sweep 
marching  for  auparaonic  outer  flows  with  little 
or  no  upstream  Interaction,  and  global  or  multiple 
peas  methods  whan  upstream  Influence,  Including 
possible  separation,  is  an  Important  feature  of 
the  flow  have  been  discussed  for  a  variety  of 
problems.  References  2  and  3  review  much  of  the 
earlier  work  on  this  subject. 

The  present  study  is  concerned  with  the  analy¬ 
sis  of  a  global  approach  to  the  Msvier-Stokss 
equations  and  a  single  sweep  PUS  application  for 
the  supersonic  laminar  flow  over  a  sharp  cone 
geonetry.  This  geometry  has  bean  examined  exten¬ 
sively  during  the  past  decade  and  in  view  of  tha 
variety  of  complex  flow  phenomena  associated  with 
Che  cone  at  Incidence,  it  has  served  as  a  prototype 
for  analysis  of  new  computational  methods.  Of 
particular  interest  here  are  procedures  that 
allow  for  the  evaluation  of  such  a  configuration 
at  relatively  large  angles  of  incidence.  In 
particular,  we  ara  concerned  with  the  accurate 
prediction  of  the  secondary  flow  vortex  patterns, 
Imbedded  shock  formation,  axial  and  secondary  flow 
pressure  interaction,  and  finally  the  surface 
pressure  and  heat  transfer. 

Tha  secondary  flow  interaction  with  shocks 


and  vortex  formation  is  non- trivial  and  therefore 
a  coupled  strongly  implicit  (CSIP)  algorithm, 
developed  previously  for  the  Incompressible  Navier- 
Stokes  equations,3  is  applied  here  In  order  to 
provide  a  "two-dimensional"  coupling  for  the  three 
velocity  components.  The  secondary  flow  presaure 
and  temperature  gradients  remain  uncoupled  in  order 
to  limit  computer  storage  requirements;  however, 
the  latter  is  evaluated  with  the  strongly  Implicit 
(SI?)  algorithn.  The  outer  shock  boundary  condi¬ 
tions  are  completely  coupled  into  the  CSIP 
algorithm. 

The  continuity  and  normal  momentum  equations 
are  related  to  the  evaluation  of  the  normal  velo¬ 
city  and  pressure,  respectively.  This  is  contrary 
to  most  conventional  Mevier-Stokes  and  PHS  proce¬ 
dures.  However,  it  is  consistent  with  the  boundary 
layer  and  thin  shock  layer  techniquea  that  have 
been  applied  In  non- interactive  analyses.  A 
deferred-corrector  is  applied  to  the  upwind  or 
boundary  layer  approximation  used  In  the  initial 
marching  sweep  for  the  axial  (O  convection  terms. 
This  K-R6  corrector  provides  a  stable  second-order 
accurate  approximation  for  tha  global  marching 
procedure. 

As  a  peripheral  aspect  of  the  present  study 
several  difference  approximations  for  the  axial 
pressure  gradient  term  pg  appearing  1m  the 
6  momentum  equation  are  considered.  This  term 
controls  the  "elliptic”  pressure  Interaction2  and 
can  produce  the  so-called  departure  Instability 
Indicative  of  the  upstream  pressure  effect.  It  is 
seen  that  with  an  appropriate  forward  difference 
approximation  to  pg,  the  departure  effect  is 
suppressed  and  a  consistent,  stable,  pressure 
Interactive  procedure,  tested  previously  for  incom¬ 
pressible  flows,  is  also  applicable  for  the  cone 
calculations.  Since  there  Is  not  strong  upstream 
Interaction  for  the  sharp  cone  geonetry  this  has 
little  effect  on  the  solutions;  however,  it  does 
provide  a  basis  for  further  consideration  with 
other  configurations. 

Solutions  are  presented  for  several  cases; 
however,  the  major  emphasis  hare  is  on  the 
Tracy  10*  half-angle  cone7  and  the  Marcillat* 

9*  half-angle  cone.  Incidence  angles  of  from 
zero  to  45*  ara  considered. 
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tamj  it  now  common 

2.  Cawralnt  Emtlo»> 

Th«  full  ataady  Navlar-Stokaa  aquae Iona  ara 
conaldarad  Cor  tha  global  relaxation  procodura. 
However,  for  tha  lnflnlta  coaa  gaonatry  conaldarad 
hara  only  a  alngla  aarchlng  eweap  la  required  aa 
upatraaa  lntaractlon  la  nagllglbla.  Tharafora 
tha  axial  dlffualon  tana  and  M  aacond-ordar 
corractor  barn  not  baan  lntroducad  In  tha  Initial 
atop.  For  othar  conflguratlona,  whara  aultlpla 
relaxation  awaapa  ara  required,  thaaa  tarna  would 
ba  lntroducad;  although,  area  than  thalr  laportanca 
may  ba  quaatloaabla. 

A  aurfaea  orlantad  coordinate  ayataa  la  coo- 
aldarad  (Fig.  1) ;  x  la  aaaaurad  along  tha  cona 
ganaratora;  y  la  noraal  to  tha  aurfaea,  and  4  la 
tha  clrcuafarantlal  or  aalauthal  angle;  4  ■  0,* 
at  tha  wind-  and  laa  planaa,  raapactlvaly. 


In  thla  coordinate  fraaa  the  governing 
ration  equation*  ara  aa  follow*: 

continuity: 
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negligible  la  the  boundary  layer  and  poeelbly  la 
the  outer  layer  (or  larger  incidence  engine. 


tint 


4*)  .} 


3x 


(Id) 


Energy-equation : 
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where  d  la  the  half-cone  angle  and  r  ■  *  alnd  + 
y  coe  d.  la  addition,  the  atate  equation  and 
viscosity- tew*  ratura  and  thernal  conductlvlty- 
tanperature  relatione  clone  the  aye ten.  The 
Sutherland  lew  la  uaed  for  u  -  p(T)  and  the 
Prandtl  nunfcer  ft  ■  u  C  /X  la  conatent. 

P 

The  eye tan  (1)  la  oon-dlnena tonal lead  with 
the  shock  layer  thlchneaa  <(x,4)  ao  that 
n  ■  y/4;  (  •  i  and  d  •  4.  The  velocities  are 
nornallaed  with  the  free  atrean  0_,  p  with  pa, 

p  with  pjlj,  T  with  U^/C^.  The  unit  laynolda 
n unbar  la  *  paDa/u„.  The  ayatan  (1)  la  than 

tranaforned  free  U,n,4)  to  (x,y,d)  coordlnataa. 
The  respective  velocities  are  (u,v,w).  It  la 
algnlfleant  that  In  the  course  of  thla  trees  for- 
aatlon 
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24  » •  4  4  to) 

(2b) 

The  final  ayatan  of  equatlona  la  quite 
formidable  and  will  be  given  in  a  separate  report. 
In  order  to  allow  for  sub sequent  eode  debugging, 
error  correction  or  future  nodi fleet Ians  In  the 
difference  equatlona  or  relaxation  procedure,  each 
tarn  la  each  equation  la  Introduced  Into  the 
coeffielenta  of  the  differenced  and  quasllln norland 
form  Independently.  Therefore  each  contribution 
la  Identifiable  and  eaally  located. 


3.  Difference  Equatlona 

The  final  tranaforned  equatlona  are  differ¬ 
enced  as  follows:  conventional  central  differ¬ 
encing  for  the  n  and  a  derivatives  In  the  axial  (() 
and  axlauthal  (a)  nonantun  and  energy  equations,  see 
figure  2  for  a  description  of  the  difference  grid. 


The  axial  convective  derivatives  J/3C  ace  upwind 
differenced,  with  a  K-K  corrector  Inserted  for 
nultlpla  sweep  applications;  l.e.. 


V"i-l  ui+l"°i 
(ut>t  -  (  A(  )FL  +  a;  d-W 


e  „  e  e 

/iayu' 
sr 


i-i 


(2FL-1) 


where  PL  ■  1  for  u^  >  0;  PL  ■  0  for  u^  <  0.  The 

superscript  (*)  denotes  solutions  fron  the  previous 
global  Iteration.  Those  are  assunsd  to  be  sere 
In  the  Initial  narchlng  weep. 

The  continuity  and  no reel  nonantun  equations, 
aa  noted  previously,  are  considered  la  the  boundary 
layer  or  thla  shock  layer  aaaaer  aa  first-order 
aquations  for  v  and  p,  respectively  and  are 
differenced  at  the  half-point  i)j+(An/2)j  with 

the  trapeaoldal  rule.  Whan  the  p^  tarn  fron  the 
n  nonantun  equation  la  uaed  la  (2a)  and  (2b), 
conventional  central  differences  are  applied  for 
the  coupled  velocity  derivatives  that  are  Intro¬ 
duced  late  the  (  and  4  nonaatna  equations. 


Although  the  pressure  will  be  uncoupled  fron  the 
(u,v,w)  CSV  algorlthn,  the  p*  and  p^  tarna  do 

couple  partially  through  the  pn  contributions  la 
(2).  This  tern  can  be  el Inina ted,  with  the  aoraal 
nonantun  aquation  (1c) ,  la  tarns  of  (u,v,w) .  In 
fact,  for  la vise id  conical  flow  4/4(  I  0  and 
therefore  we  would  expect  the  role  of  p,  to  be 
relatively  ndaor.  The  other  (-derivatives  are  not 


The  final  ays ten  of  difference  equations  ate 
quasl-llaaarlsed  to  provide  the  fanlllar  second- 
ordar  approximation 


where  tha  superscript  (*)  here  denotes  the  previous 
Iteration  at  the  axial  (g  location.  Strictly 
speaking  If  (.,  values  are  uaed  this  approximation 


rasa  In*  second-order  and  la  Chan  non-ltaraelva 


algorltha  ta  atructurad  to  evaluate  all  k  values 
before  p roc aa  ding  fron  J  to  J+l. 


* 


4.  Hunerical  Alaorlthn  and  Solution  Procedure 

Altar  differencing,  aubatltutloo  of  pn  and 
quaai-llnearixaeioa ,  tha  final  •yatas  to  b« 
aval ua tad  la  of  tha  following  torn  for  (u, v,w) 
at  C  -  qs 


€.4  aonantun: 

E.  .  +  W.  U.  ,  .  ♦  H,  0.  ... 

tuJ,k  J+1,k  *ttJ,k  3_1,k  tuj,k  J,k*1 


+  S.  0.  .  ,  +  P,  V.  .  +  I.  W. . 

luJ,k  J,k"1  lttJ,k  J,k  WJ,k  2+1,k 


This  la  a  (2x2)+l  aatrix  systea  and  not  a 
such  largar  (3x3)  systea.  Tha  algorlthaa  (4a)  and 
(4b)  ara  uaad  in  (3a,b)  to  allnlnata  k> 

Vl.k*  Di.k-1*  H3,k-1*  Th*  con“t  3tl*  ktl 

ara  traatad  explicitly  In  obtaining  tha  ra cur a ion 
ralationahipa  for  tha  coafflelanta  in  (4a),  (4b). 
Saa  rafaranca  5  for  a  aatrix  description  of  tha 
CSIP  algorltha.  Thaaa  ara  obtalnad  by  dlraetly 
relating  (3a,b)  to  (4a,b).  Tha  final  recursion 
ralationahipa  ara  given  In  Appendix  A.  Tha 
boundary  conditions  for  u,  v,  w,  to  ba  discussed 
in  tha  following  aectlon,  ara  coupled  Illicitly 
into  tha  CSIP  algorltha. 


♦  W 

+  P 
+  E 


J.k 


J.k 


tv 


J.k 


3-1, k 


J.k 


J+l,k 


+ 

+ 

+ 


"^.k 

Wl'j.k  Vw 


+ 


+ 


VJfkWJ,k-l 

P»vJfkVJ.k 
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Tha  energy  aquation  which  la  of  tha  fora  (3a) 
la  solved  Independently  with  tha  SIP  algorltha  (4a). 
This  la  obtalnad  by  replacing  u  with  T  and  setting 
v,w  -  0  In  (3a)  and  (4a)  and  Inserting  the 
appropriate  coefficients  froa  tha  transforaed  (la). 
The  pressure  Is  obtalnad  froa  the  tranaforaed 
noraal  aoaantua  aquation  (lc)  by  Integrating  froa 
tha  shock  to  tha  surface  with  tha  trapesoldal  rule. 


where  1  »  2,  3  for  C,4  aquations  respectively, 
continuity: 

VJ.k  "  WlvJ  k  Yl.k  +  Pluj  k  uJ,k 
+  V  .  DJ-l,k  +  Plw  .  Hj,k 

Jt* 

1Wj,k  i“1»k_1  lj,k 

.0  <»> 

The  coupled  strongly  lapllelt  (CSIP)  algorltha 
Is  of  tha  tort? 

°J.k  "  *llj  k  Vl.k  +  *12^  DJ,k41 
+  *11,  k  Vl.k  +  *12,  .  WJ,k+l 

*Vw‘V  . 


"j.k  "  *21^  Vl.k  +  *22j  k  DJ,k+l 
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(4b) 


sad  tha  conelauley  aquation  (3b)  where  tha  comer 
polata  ara  treated  lapUcitly.  Tha 


Tha  procadura  can  ba  suaaarlsed  as  follows: 

(1)  (u.v.w)  ara  solved  with  the  CSIP  algorltha 
and  coupled  boundary  conditions  with  p,  T,  p 
prescribed  (pp  has  bean  substituted  froa  aoaantua 
Into  tha  C,4  aoaantua  aquations  and  is  thereby 
lapllelt);  (2)  With  the  new  shock  values  for  p 
and  T,  T  la  updated  with  tha  SIP  algorltha  and 
p  froa  tha  aornal  aoaantua  aquation.  Tha  viscosity 
and  thexaal  conductivity  coefficients  and  p  ara 
also  updated;  (3)  With  tha  new  p,  p,  T,  u,  k 
values,  return  ta  step  (1)  and  continue  until  a 
convergence  criterion  is  satisfied.  A  flow  chart 
for  tha  entire  procadura  Is  given  In  Appendix  1. 

This  describes  caaplately  tha  evaluation  at 
each  location  (  «  (j  for  a  single  sweep  nsrchlng 
procedure.  If  additional  sweeps  are  necessary, 
l.e. ,  whan  pressure  Interaction  and  upstrean 
Influence  are  laportant,  the  process  is  repeated 
fron  the  Initial  (  *  C«  location  and  the  EX 
correctors,  axial  diffusion,  forward  p,  difference 
tamo  dependant  on  the  previous  pass  are  updated. 
Multiple  pass  calculations  have  only  been  consi¬ 
dered  to  a  United  extent  and  these  will  not  be 
discussed  further  here. 


5.  Boundary  and  Initial  Conditions 

Initial  conditions  arc  obtained  directly  fron 
the  governing  difference  aquations  by  applying 
the  siasrlcal  algorltha  with  i/H  »  0;  l.e., 
a  conical  flow  approxlnatloo.  Invlscld  or  step 
profiles  with  no-sllp  surface  conditions  ate 
■Miami  to  start  this  calculation.  This  procedure 
works  well  for  Incidence  angles  up  to  10*.  Ibi 
a  >  10*,  the  10*  results  are  used  as  Initial 
conditions  and  the  exact  a  boundary  conditions 
are  loosed  at  the  first  narehlng  step.  The 
effect  of  the  Initial  conditions  rapidly  vanishes 
as  the  shock  layer  Is  extreasly  mall  near  the  tip 
of  the  erne.  It  should  be  noted  that  while  the 
}/>(  ■  0  solutions  are  reasonable  approxlnatlona 
for  the  viscous  cone  flew,  the  exact  solutions  are 
noticeably  different  both  In  the  viscous  boundary 
layer  and  la  tha  outer  vertical  region  near  the 


1m  plane.  This  la  net  surprising  as  the  vlacoua 
behavior  does  not  aeala  conically,  i.a. ,  aa  i(x,t) 


Ths  boundary  condltlona  at  tha  surface  n  •  0 
ara  u  ■  v  •  w  ■  0  and  T  •  1^.  At  tha  outer  ahock 
vava  n  ■  1,  the  lank 1 na-Hugonlot  condltlona 
apply: 

aaaa  conaarvatlon: 

(Uj-  pu)«e  i  (Wj-  pw)d^  -  (vx  -  pv)  -  0  (5a) 


tangential  noaantun: 

(Wj-w)  +i  (Vx-V)  -  0  (5b) 


(uj-u)  +  «?  (Vj^- V)  -  0  (5c) 


noraal  noaantun: 


p  - - **■*■■■=—  (5d) 

T-1  +  -J- 

\ 

abate  aubecrlpt  (1)  denotea  free  a trees  conponenta 
at  the  shock  wave  and  H,,  la  the  free  a  trees  noraal 
Msch  nuaber  which  la  a  x function  of  i^. 

The  four  equatlona  (5a-d)  coupled  with  tha  con¬ 
tinuity  equation  (la)  evaluated  at  the  half-point 

give  five  equatlona  for  the  five 

uaknowna,  u,  v,  w,  p,  The  direction  coelae 
la  not  coupled  but  la  treated  Iteratively  with 
p  and  T  which  ara  only  fmctlooa  of  1^.  fron  the 
Hawhlna  Bugoniot  conditions.  The  coupled  u,  v,  w 
botmdary  condltlona  are  Incorporated  lapllcltly 
Into  the  CSIP  algorithm.  The  ayatea  (5)  la  quasi- 
llnearised  la  the  aaaa  wanner  aa  that  for  the 
Interior  equatlona.  The  coablnatlou  of  tha 
lanklne-Hugoniot  condltlona  and  the  coupled  con¬ 
tinuity  equation  enaure  that  global  aaaa  conser¬ 
vation  la  smlnta load.  The  final  torn  of  the 
outer  boundary  conditions  lncludea  (5b) ,  which  la 
linear  alnce  la  treated  explicitly,  and 
(5e,c,d)  are  cosblnad  to  give 

j  [1-  (UjO+  VjV+  WjW)  ]  ( (r-1) 

♦  (2-  ) (UjU+  v^v+  WjW)  +  (u2  +  v2  +  w2)  ] 

•  II-  2(u1u+v1v+WjW)  +  (u2+v2+w2)1  (5e) 


This  equation  la  also  quasl-llnaarlted.  In  this 
fern  the  strean  ftmctlon  solution,  appearing  la 
the  original  ayatea,  has  been  eliminated. 

This  CSX?  and  SI?  procadures  are  convergent 
to  10-1°  for  all  Incidence  angles.  Tha  pressure 
equation  la  converged  to  10 “7  for  a  _<  32*.  Per 
a  >  32*  the  pressure  solution  exhibits  a  slow 
divergence  that  la  apparently  associated  with  the 
explicit  or  tacoupled  character  of  the  nomal 
noaantun  aquation.  If  tha  praasure  field  la 
tmdarrelaxsd  by  Including  a  tins- like  tern  la  the 


noraal  aoasntua  and  energy  equations,  this  diver¬ 
gence  la  suppressed  and  tha  overall  systea  con¬ 
verges;  however,  the  uaderrelaxatloa  la  very  expen¬ 
sive  coaputatlonally  and  It  la  evident  that 
coupling  the  pressure  directly  would  be  aore 
desirable. 

The  solutions  for  a  >  30*  depict  fairly  strong 
labedded  ahock  waves  and  significant  vortex  lift¬ 
off;  therefore,  the  secondary  flow  pressure  inter¬ 
action  will  be  quite  strong  and  an  Implicit  txeat- 
aent  of  this  effect  nay  be  required.  He  recall 
that  the  p.  contribution  la  (2b)  Is  already  Intro¬ 
duced  Illicitly;  however,  the  p^  and  if  effects 
are  considered  Iteratively. 

finally,  sywetry  conditions  u^  •  v^  •  ■ 

p^  •  0  and  w^  ■  0  are  imposed  at  the  lee  and 

wind  planes.  These  conditions  are  coupled  lapllclt¬ 
ly  Into  the  CSX?  and  SI?  algorlthaa  by  appropriately 
aedlfylag  the  coefficients  of  the  difference 
aquations  at  these  locations;  e.g. ,  uk+1  •  ufc_1 

at  k  ■  1,  therefore  S^  •  0  la  (3a,b)  and  Htu^k 
la  nodi  fled  to  Include  the  0^^  coefficient  so 
that  R|uj|  "  fluji  ^tuj ^ •  etc. 

6.  Pressure  Interaction  and  Departure 

Iter  any  aweep  of  a  one -step  or  global  relaxa¬ 
tion  procedure  the  appropriate  treatment  of  the 
pressure  Interaction  term  Pg  la  the  axial  nomen  tun 
equatlona  (lb)  la  eaaentlal.  Sven  for  supersonic 
free  streams  am  "elliptic"  pressure  Interaction 
la  transmitted  upstream  through  the  term  rinn 
regions  of  subsonic  flow,  as  la  the  wall  beivdary 
layer,  are  present.  If  Pg  la  prescribed  tha  equa¬ 
tlona  axe  truly  parabolic  and  the  upstream  effect 
la  lost.  However,  for  flowe  where  upstream 
effects  axe  aip>lf leant,  if  the  p(  term  la  backward 
differenced  In  any  fora  the  Interaction  la  included 
i'r  '  — 

In  tha  ayatea.  '  This  la  manifested  la  the 
appearance  of  departure  curves  or  Instabilities 
whenever  Ax  >  (Ax)  which  la  of  the  order  of 

“  “  -I/? 

tha  thickness  of  the  subsonic  layer  or  0(Be  ) 
for  tha  vlacoua  eoea  flow.*  Therefore,  In  the 
present  application,  with  He  »  1,  reasonably 
accurate  eolations  are  possible  even  for  Ax  ■  (Ax),,,. 

This  procedure  la  however  inconsistent  and 
for  flows  with  siplficamt  pressure  Interaction 
nay  be  Inadequate.  A  global  relaxation  procedure 
la  required  If  tha  inconsistency  la  to  be  elimi¬ 
nated.  Me  rover,  la  the  global  procedure  It  la 
necessary  that  each  marching  sweep  be  departure 
free  and  yet  allow  for  an  Implicit  pressure  inter¬ 
action  la  order  to  circumvent  the  separation 
point  singularity  that  arises  with  sm  explicit 
pressure  representation.  In  an  earlier  incompres¬ 
sible  TVS  thin  layer  study,***  it  was  shown  that 
forward  pe  differencing  la  departure  free  and 
stable  globally.  Central  p,  la  nestable  globally 
and  backward  pc  la  not  departure  free;  It  was 
shown  la  references  2  and  •  that  with  forward  P( 
differencing,  It  wee  possible  to  aareh  through 
separation  regions  departure  free  la  each  global 
relaxation  sweep.  This  procedure  was  also  tasted 
for  tha  cone  calculations  and  the  earlier  conclu¬ 
sions  obtained  both  analytically  and  numerically 
ware  verified.  Although  this  affect  Is  relatively 


a 


% 


tai  important  for  cob*  flows  whoro  solutions  with 
As  >  (ii)^  or*  KMfUbU  except  wry  clot* 

to  tho  tip,  for  mn  complex  configurations  this 
procedure  is  considered  l^ortant. 

For  tho  coos  calculations ,  da  tall ad  tasts  for 
tho  value  of  (Ax)..-  as  a  function  of  la  haws 
boon  considered.  Toe  resulting  stability  curva 
for  o  “  21*  la  shown  in  figure  3.  This  clearly 
confine  tho  oarliar  analysis  sad  dsnonstratas  that 
(Ax)  ,_/y„  ,  dacroasos  aa  la  increases.  Tha  exact 
ala  wi  .  12 

dapondanco  Is  not  proportional  to  la  aa  In  tha 
Incompressible  case.  Tho  use  of  a  l-l  second-order 
pressure  correction  for  a  forward-differenced  px 
has  not  bean  addressed.  Finally,  we  recall  that 
the  contribution  to  p_  la  treated  inpllcitly  and 
with  a  backward  or  merchug  nods  for  v,;  therefore, 
the  entire  global  Interaction  Is  contained  in  pj. 


JhU 

^•"Ul 


Fig.  3.  Departure  Condition  with  P  Backward 
Differencing.  ' 


f taction  of  a  to  reach  the  sans  C  location  (•  2  ft) 
for  Tracy  conditions.  Free  the  figure  It  la  seen 
that  running  tins  Increases  with  a  and  due  to 
imdarrelaxatloa  required  for  a  >  32*,  the  tines 
lacrosse  Measurably.  Double  precision  requires 
approximately  four  tines  that  of  single  precision 
which  runs  about  four  tines  slower  than  a  CDC  7600. 
Therefore  we  estlaete  conservatively  that  the 
naxlnun  CPU  requirement  on  a  7600  would  be  less 
than  25Z  of  that  of  the  AMDAHL  co^uter  If  similar 
prograaning  techniques  were  used.  Non -unifora 
grids,  in  particular  for  the  surface  no real  (n) 
direction,  night  reduce  the  overall  cosputar  tins; 
however,  for  large  incidence  with  labedded  shocks 
and  lifting  vortex  notion  the  optimal  grid 
distribution  Is  not  obvious. 


Fig.  A.  CPU  Tims  for  Initial  (jr  -  0)  and  Marching 
Solution.  4 

Solutions  have  boon  obtained  for  incidence 
angles  (a)  up  to  65*  on  a  10*  half-angle  (•)  cone, 
or  a  :  6.5  8.  The  Interior  shocks  are  captured, 
secondary  flow  separation  regions  are  determined 
and  vortex  patterns  are  seen  to  lift  far  off  the 
surfacs  for  large  o. 


7.  Basulta 

Solutions  have  been  obtained  for  two  flow 
conditions  for  which  experimental  data  is  avail¬ 
able'  and  for  which  numerical  results  have  bean 
obtained  by  merged  layer,  “  and  boundary 

region^  approximations .  All  tbs  calculations 
were  mads  with  the  C8IP  algorithm  for  (u,v,w),  the 
SIP  algorithm  for  T  and  the  trapesoldal  rule  for  p. 
Converged  solutions  to  a  tolerance  of  U*“  for 
successive  Iterations  using  double  precision  corn- 
p lletlon  on  aa  AWIAHIi  370-computer  were  obtained. 
Uniform  gride  were  specified  la  u  and  6.  Calcu¬ 
lations  ware  nai  with  both  relatively  coarse  grids, 
e.g.,  10  x  21,  19  x  31  and  a  finer  51  x  55  grid. 

For  the  finest  grid  calculation  tints  wars 
0.0008  sec/CSIP  Iteration/grid  point.  7-15  CSIP 
Iterations  era  required  to  converge  the  (u,v,w,4.) 
system.  The  CSIP  procedure  Is  called  about 
30-60  times  to  converge  the  pressure  field.  FOr 
s  32*  it  takas  about  30  minutes  for  the  (»/»?■  0) 
equations  to  converge  ad  another  1  minute  for 
any  additional  marching  step. 


Selected  results  are  presented  here  primarily 
for  the  experimental  conditions  of  Tracy7  and 
Mar  cl  list.’  Sons  of  tbs  solutions  are  for  very 
large  incidence  angles,  where  experimental  data 
are  ixia  vail  able.  A  51  x  55  uniform  grid  Is  the 
finest  mesh  considered  to  date  and  therefore,  we 
Inject  a  note  of  caution  for  o  >  32*.  In  addition, 
as  discussed  previously,  underrelaxation  Is 
required  for  these  cases  and  further  study  of 
accuracy  and  convergence  properties  is  required. 

The  heat  transfer  at  the  wind  and  lee  planes 
for  moderate  to  large  angles  of  incidence  is 
depleted  in  figures  5e,b.  The  agreement  with  the 
data  Is  quits  good  as  It  is  for  the  surface  pres¬ 
sure  distributions  which  are  given  in  figures 
6a,b,e.  The  effect  of  the  labedded  secondary  flow 
shock  wavs  Is  apparent  only  for  the  largest  inci¬ 
dence  angles. 

These  affects  ean  be  seen  somewhat  more 

clearly  in  the  a  ■  26*  pitot  pressure  profiles 
and  distributions  of  figures  7a,b,c.  The  signi¬ 
ficant  change  la  p0  across  the  shock  layer  Is 
apparent.  The  stagnation  pressure  is  affected  by 
the  Inclination  of  tha  outer  shock,  the  growth  of 
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Fig.  5a.  Syunetry  Plana  Haat  Tran* far  - 

Moderate  Incldanca.  Fig.  6a.  Surfaca  Praaaura  at  Moderate  Incidence. 


Fig.  Sb.  Syanetry  Plana  Haat  Traaafar  - 
Larga  Incldanca. 


tha  boundary  layer,  aacondary  flew  aaparatlon  and 
the  appaaranca  of  lnbaddad  Chock  wavaa.  Bear  tha 
w lad  plana,  ehaaa  propartlaa  an  aaaant tally 
Invariant  aa  la  p0-  However,  aa  la  aaan  la  tha 
following  flguraa,  tha  aacondary  flew  pattarna 
fron  tha  wind  plana  an  axtranaly  complex  and  In- 
cluda  all  of  tha  phanonana  alludad  to  pravloualy. 

Tha  aurfaca  atraaallaa  Inclination  and  flow 
pattarna  (figuna  I,  9  and  Tabla  2)  raflact  tha 
occurranca  of  aacondary  flow  aaparatlon  and  for  tha 
largaat  Incldanca  angina  a  double  aaparatlon 
affact  la  pradlctad.  Thla  phanonana  haa  boon  aaan 
azparlaantally  In  turbulent  f  Iowa  and  by  Mardllat9 
for  lanlaar  flow*,  but  haa  not  boon  nportad  In 
carllar  nuaarlcal  atudlaa.  In  view  of  tha  rala- 
tlvaly  coaraa  51  x  55  uniform  grid  wa  raaarva 
conaant  on  tha  accuracy  of  the  praaant  raaulta; 
howavar.  tha  vary  axiatanca  of  tha  convarged  doubla 
aaparatlon  aolutlona  obtained  hara  la  lntaraatlng. 
Tha  primary  aaparatlon  point  la  ralatlvaly 
lnaanaltiva  to  Incldanca  for  larga  a  and  la  ahown 
In  figure  10.  Aa  la  all  aarllar  calculatlona, 
for  tha  praaant  analyala  aaparatlon  flrat  occura 
for  o  :  7*  -  8*  and  tha  aaparatad  raglona  found 
by  Merclllat9  for  a  >  6*  ara  not  racovarad  hara. 

Tha  location  of  tha  prlaary  vortax  cantor  la 
depleted  la  figure  11.  The  vortax  aovaa  away  fron 
tha  aurfaca  for  lacraaalag  «,  but  appeara  to 
aaynptoto  to  a  fined  location  at  a  given  value  of 


Fig.  6e.  Surface  Praaaura  at  Vary  Large  Incldanca. 
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JODIS  (1968) 
(lnvlscld) 


(55  x  51  Grid) 


JODIS  (1968) 
a-U«  (1“Tl‘cl<1) 

M_  - 10  LID  6  IUBIN  (1973) 
6-10* 

present  method 
(55  x  51  Grid) 


>  22.5 

65 

67.5 

90 

112.5 

135 

157.5  1 

>  1.732 

3.286 

6.583 

5.517 

5.851 

5.385 

3.502 

l  1.726 

3.291 

6.601 

5.533 

5.870 

5.390 

3.671 

i  1.722 

3.288 

6.598 

5.531 

5.872 

5.392 

3.697 

i  2.667 

5.333 

7.017 

9.217 

9.300 

10.300 

7.239 

I  2.830 

5.701 

7.889 

9.650 

10.619 

11.090 

8.050 

I  2.896 

5.852 

8.113 

9.938 

10.908 

11.686 

8.696 

TABU  2.  SURFACE  STREAMLINE  INCLINATION,  tan  (v  /u  )| 
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Fig.  9.  Surface  StrtuUM  Isttim,  »■  12*.  26*. 
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fig.  10.  Prlnar;  Separation  Folnt  Location. 


Tig.  11.  Trlnery  Vortex  Cantor  Location. 

t\  relatival;  cloaa  to  tha  shock  tor  o  *  30*.  Tha 
"inviacid"  aolutlon  In  the  outer  portion  of  tha 
ahock  layer,  Including  tha  vortax  location  la 
sp proxiastel;  lnda pendant  of  C  and  a  function 
only  of  i)  and  6,  l.e.,  conical  flow  behavior  la 
obaarvnd  there. 

Tha  boundary  and  labedded  ahock  patMiu 
are  aeon  la  flgurea  12a, b,c  for  aavnral  Incidence 
an  glee.  Tha  labedded  ahocka  Incraaaa  In  atmngth 
and  tha  outer  ahock  ehape  begins  to  deviate 
alpilflcantly  fron  tha  nonr  circular  pattern  at 
•aallar  Incidence  angle  a.  Dm  flow  la  super- 
aonlc  over  aoat  of  tha  mage  of  6  for  a  >  12*. 

finally,  for  a  ■  65*,  tha  variation  of  the 
ayanatry  plana  surface  pressure  and  a  qualitative 
picture  of  tha  eecoedary  flow  la  glean  in 
flgurea  13a,b.  Tha  strong  labedded  ahocka,  tha 
double  vortex  pattern  and  tha  apple  shape  of  tha 
outer  ahock  boundary  am  of  particular  note. 

These  solutions  am  converged  to  the  tolerances 
described  earlier;  however,  the  grid  any  be  too 
coarse  to  reach  a  final  conclusion  on  the  accuracy 
of  this  pic turn. 


Fig.  12a.  Shock  Pattarns  -  a  -  12*. 


Fig.  13a.  Surface  Praaaura  on  SyMetry  Plana. 


♦  svemo  mnn 

Fig.  13b.  Flow  Pattern  for  a  -  43*. 


Fig.  12b.  Shock  Patterns  -  a  -  24*. 


In  eaeiry,  the  CSIP  procedure  deecrlbed 
here  captures  ehocke,  predicta  vortex  lift-off, 
Multiple  separation  regions  and  conplex  Shock 
patterns.  The  possible  suppression  of  departure 
effects  with  global  iteration  and  forward  pg 
differencing  and  the  adequacy  of  single  sweep 
backward  differencing  and  the  adequacy  of  single 
sweep  backward  differencing  for  conical  flown 
with  snail  subsonic  viscous  layers  has  been 
denonstrated.  The  accuracy  of  the  present  solu¬ 
tions  for  vary  large  incidence  ranslns  a  question 
for  further  study. 
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